Twisting and buckling: a new undulation mechanism for artificial swimmers 
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We present an artificial swimmer consisting in a long cylinder of ferrogel which is polarized trans- 
versely and in opposite directions at each extremity. When it is placed on a water film and submitted 
to a transverse oscillating magnetic field, this artificial worm undulates and swims. Whereas sym- 
metry breaking is due to the field gradient, the undulations of the worm result from a torsional 
buckling instability as the polarized ends tend to align with the applied magnetic field. The critical 
magnetic field above which buckling and subsequent swimming is observed may be predicted using 
elasticity equations including the effect of a magnetic torque. As the length of the worm is varied, 
several undulation modes are observed which are in good agreement with the bending modes of an 
elastic rod with free ends. 

PACS numbers: 47.63.Gd; 47.63.mf; 46.32.+x; 46.40.-f 
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The design of robots that mimic animal locomotion has 
gained in interest during the last two decades, both for 
their medical applications as well as for exploration tasks 
or search and rescue missions. Among the various loco- 
motion strategies of animals, the undulation locomotion 
is of particular interest for its simplicity and robustness, 
even on rough terrain. The locomotion of snakes, eels or 
worms have been particularly studied from fundamental 
as well as biomimetic points of views [l|-[3| and robots 
inspired by their displacement methods have been pro- 
posed (1,|5|. To wirelessly power and control such robots, 
smart materials have been developed, often actuated by 
magnetic fields. The undulation of artificial swimmers 
are generally obtained by oscillating a magnetic head at- 
tached to a flexible tail j6|-|l|. Beyond the technological 
challenge, such swimmers also raise a more fundamen- 
tal interest on the description of the non-reversible flow 
around the body |9l-[T2j. In this paper, we present a new 
and non-trivial undulation mechanism that is based on 
the torsional buckling instability of a worm-like swimmer 
under an oscillating magnetic field. 

Our swimmer consists in a flexible cylinder of 
PVA (poly- (vinyl- alcohol), Sigma- Aldrich) based ferro- 
gel, reticulated with glutaraldehyde (Sigma- Aldrich), in 
which 10% in weight of ferromagnetic particles (black 
iron oxide), with micrometric size, are added. The parti- 
cles are embedded in the gel, which has a Young modulus 
E = 1.36 kPa (measured by indentation on a TA-XT2 
texturometer). Before reticulation, the gel is poured into 
a cylindrical mold, with diameter 2R = 1 mm and length 
10 < L < 50 mm. Once the gelation is achieved, the 
gel is taken out of the mold and rinsed with water to 
remove free PVA chains. In order to provide perma- 
nent magnetic properties to the swimmer, a magnet is 



then placed at each of its extremities, with opposite po- 
larization directions. Over a distance a = 5 mm from 
each extremity, the embedded particles acquire a perma- 
nent magnetization density M, oriented perpendicularly 
to the main axis of the cylinder, with opposite direction, 
see Fig. [TJa). The magnetization density was measured 
on the magnetization curve according to Foner's method 
and equals 1.5 G. The magnetic worm-like swim- 
mer is then deposited at the surface of water in a Petri 
dish, set between two coils in Helmholtz configuration 
(Figlljb)). When we supplied with AC, the coils gener- 
ate a vertical oscillating magnetic field, with adjustable 
magnitude B and frequency /. 
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FIG. 1: (a) Sketch of the magnetic worm-like swimmer: the 
colored zones represent the polarized regions, (b) Experimen- 
tal set-up: a flexible magnetic worm is placed at the surface of 
water, between two coils in Helmholtz configuration, generat- 
ing a vertical magnetic field B. (c) Top view of the swimmer 
under an oscillating magnetic field (B = 50 G, / = 1 Hz). The 
images are thresholded for better visualization. The time step 
between successive images is 0.33 s, and the length L of the 
swimmer is 35 mm. 

Submitted to a constant magnetic field {i.e. for / = 
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Hz), the worm drifts slowly from the center to the edge 
of the Petri dish, without deformation 0|. The dis- 
placement velocity V is very small 0.2 mm.s~^) and 
proportional to the magnitude of the magnetic field 5, 
Fig. [21 The motion is due to the existence of a gradient 
of the horizontal component of the magnetic field in the 
radial direction, which increases by 30 % from the center 
to the edge of the dish. A magnetic force, proportional to 
the field gradient, and thus to the field magnitude, acts 
on all the ferromagnetic particles of the worm (not only 
its extremities) and pushes it away from the center of the 
Petri dish. 

Under an oscillating magnetic field, the worm undu- 
lates (Fig. HJc)) and swims toward the edge of the dish 
with a much higher velocity, see This only occurs 

above a critical value of the applied magnetic field. Be 
(Fig. [2]). Below this value, the swimmer drifts slowly 
without undulating, just as when subjected to a constant 
magnetic field. 
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FIG. 2: Swimming velocity V of the worm as a function of 
the magnitude B of the applied magnetic field, for a constant 
field (squares) and for an oscillating field at / = 1 Hz (circles). 
The length of the swimmer is L = 15.6 mm. The critical field, 
Be, is here 40 G. Inset: maximum curvature k of the swimmer 
as a function of the magnetic field, for f — 1 Hz. 

In order to investigate the origin of the undulation, 
white dots were drawn along the worm body to serve as 
labels. The dots allow to highlight the full 3D deforma- 
tion of the swimmer and reveal that the worm actually 
twists around its main axis when it undulates [14]. The 
magnetic moments of the extremities tend to align with 
the applied vertical magnetic field, and since they are po- 
larized in opposite direction, this results in the twist of 
the ferrogel cylinder. Above a certain threshold, twisting 
of an elastic rod is known to lead to its buckling through 
an elastic instability, this is known as the Greenhill prob- 
lem : a twisted elastic rod buckles into a 3D helix. As 
the magnetic swimmer is here constrained in the plane 
of the free surface of the liquid by surface tension, the 
deformed shape of the worm is a confined helix, a quasi- 
sinusoid [J^] . The alternative twisting of the worm under 



the oscillating magnetic field would therefore be at the 
origin of the bending wave. 

To test this hypothesis, the critical field Be was mea- 
sured for worms of different lengths (Fig. [3|). For an 
elastic rod with Young modulus E and twist modulus C, 
the total elastic energy reads: 
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with / = 7ri?^/4 the area moment of inertia and s a curvi- 
linear abscissa along the rod {w]. The first term corre- 
sponds to the bending energy, being the local bending 
angle, and the second term is the torsion energy, with 
(f) the local twisting angle. Both terms are proportional 
to EI/L. The rod is submitted to a twisting torque due 
to the interaction of the magnetic moments of the ex- 
tremities with the applied magnetic field, whose modu- 
lus reads T = \\7r R'^aJ^ A ^||, with irR'^a the volume of 
the polarized region, and M the magnetization density. 
Comparing the elastic and magnetic energies shows that 
the critical twisting torque above which buckling occurs 
should scale as EI/L and thus the critical magnetic field 
Be as ER^ / {aML). A more detailed description of the 
problem, based on the elasticity equations coupled with a 
magnetic torque, allows to obtain the missing pre-factor 
in the above expression. The critical field is then found 
as [HJ: 
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This model is compared to the experimental data of Be 
as a function of the swimmer's length, see Fig.[3l A good 
qualitative and quantitative agreement is found, support- 
ing that the twisting of the worm due to the interaction of 
its polarized extremities with the applied magnetic field 
is indeed at the origin of the undulation of the swimmer. 
The scaling of the swimming threshold. Be ~ R^/aL, 
is favorable to miniaturization since it is globally inde- 
pendent on the swimmer's dimensions. The propagative 
nature of this bending wave, allowing the worm to swim, 
originates from the asymmetry of the twisting torque on 
both extremities, resulting from the radial gradient of the 
applied field. This gradient is here inherent to the exper- 
imental set-up. Another way of breaking the head/tail 
symmetry is to unbalance the magnetizations of the ex- 
tremities (see p3|). 

Above the critical field Be^ the swimming velocity of 
the worm increases with the applied field, with a slope 
twenty times larger than below Be (Fig. [2]). This increase 
is associated with a rise in the body curvature with B, 
Fig. [2] (insert). As the total elastic energy EI/L) is 
inversely proportional to its length, this effect is all the 
more pronounced as the length of the swimmer is larger 
(Fig. El insert). The increase in the slope /SV//S.II with 
L is in fact due to the sum of two effects: the increase 
in the fiexibility of the rod with L, and the fact that 
the gradient of the magnetic field and thus the difference 
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FIG. 3: Critical field Be as a function of the length L of the 
swimmer (/ = 1 Hz). Dash line: fit of the data, of the form 
Be = a/L {a = 649 G.mm); solid line: model, Eq. (2] using 
the experimental values of parameters (no adjustable param- 
eter): E = 1.36 kPa, R = 0.5 mm, a = 5 mm, M = 1.5 G. 
Insert: variation of velocity beyond the threshold with the ap- 
plied magnetic field as a function of the length of the swimmer 
(dotted line: linear fit). 



of magnetic torque applied at both extremities (at the 
origin of the symmetry breaking) also increase with the 
length of the swimmer. 
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FIG. 4: Curvature of the worm as a function of its length, for 
B — 50 G and / = 1 Hz. The inserts are the superpositions of 
the successive shapes of the worm over one period, reveaUng 
the successively obtained mode shapes. 

Interestingly, as the length of the swimmer is increased 
(for a fixed value of B and /), different undulation modes 
are observed and characterized by n, the number of antin- 
odes, Fig. HI These successive modes of deformation are 
associated with an increasing curvature for increasing 
L. Note that modes 2 and 3 cannot be explained by 
a static twisting-buckling instability as they would re- 
quire a twisting angle larger than tt which, in practice, is 
impossible here: if the magnetic moments of the extremi- 
ties get aligned with the applied field (which happens for 



a torsion angle of tt), the magnetic torque cancels and 
no further rotation can happen. Another type of expla- 
nation has therefore to be sought through the dynamic 
behavior of the swimmer. 

The bending amplitude A and swimming velocity V 
are measured as a function of the frequency of the applied 
magnetic field, for fixed values of L and 5, Fig.O Both A 
and V increase approximately linearly with / up to 5 Hz, 
where they reach a maximum, and then decrease. The 
curve of the amplitude displays a pronounced peak which 
resembles a resonance curve. The peak frequency may be 
compared with the free vibration frequencies of an elas- 
tic rod. Two main modes of vibrations could be relevant 
here: bending modes and torsion modes (or coupling be- 
tween both). For bending modes, the boundary condi- 
tions for a rod with free extremities in mode 2 (the mode 
observed at all frequencies in Fig. [5]) impose kL = 7.85, 
with k the wave number. Using the dispersion relation 
uj = (3k'^ {(3 = ^EI/{piTR^), with p = 1050 kg.m-^ the 
mass per unit volume of the rod), the eigenfrequency is 
found as IITI 
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Using the previously mentioned values of the param- 
eters (and L = 29 mm), /n=2 is here equal to 3.44 Hz, 
close to the experimental peak frequency at 5 Hz. The 
eigenfrequency of a torsion rod, for n = 2, is of the or- 
der of 230 Hz: the torsion modes alone are therefore 
not relevant to account for the dynamics of the mag- 
netic swimmer. However, coupling between the bending 
modes (which are predominant) and the torsional modes 
is needed to explain that the peak frequency is higher 
than the eigenfrequency of pure bending. 
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FIG. 5: Bending amplitude A (a) and velocity V (b) of the 
swimmer (the / = Hz drift velocity, 0.23 mm.s~^, has been 
substracted) as a function of the frequency / of the oscillating 
magnetic field (B = 50 G, L = 29 mm). Solid line: model 
(Eq. |4| with the sole mode 2; dash line: model with mode 2 
up to 10 Hz and mode 3 above 15Hz. 

The evolution of the swimming velocity with the ap- 
plied frequency may now be explained thanks to a classi- 
cal kinematic model. The cutting frequency between in- 
ertial and viscous regime, fc = 2ri/ pur'^ (with r] = 10~^ 
Pa.s~^ the viscosity of water), is here of the order of 2.6 
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Hz: over the frequency range studied here, we are there- 
fore in a mixed regime where viscous and inertial effects 
coexist. As the progression of the worm happens at rel- 
atively low Reynolds numbers (from 0.3 to 2 and from 
1 to 20, respectively for the longitudinal and transverse 
motion), we use the predictions of the swimming veloc- 
ity in the Re ^ limit although the bending motion of 
the swimmer is mainly dominated by inertia. Balancing 
the longitudinal forces acting on a sinusoidal elongated 
worm, submitted to purely viscous friction, allows to ex- 
press the progression velocity as [18|-[2Q| 

^^. AVi,^ ,4, 

^ [2^ A^k^{c- 1/2) j 

with A, u = 27r/, and /c, respectively the amplitude, 
pulsation and wavenumber of the propagative wave, and 
c the ratio of the transverse to longitudinal friction co- 
efficients (c = 2 [U, HH). In the above equation, k is 
obtained from the wavelength measured on the superpo- 
sition of the successive shapes of the swimmer (as shown 
Fig. Hj). The main mode of deformation is mode 2 for all 
frequencies up to 10 Hz, but a combination of modes 2 
and 3 is observed above 15 Hz, leading to two possible 
values of the wavenumber in Eq. [H at high frequency. 
The amplitude A to be considered in Eq. [4] is that of 
the propagative component of the bending wave animat- 
ing the worm: the wave is indeed mainly stationary in 



our case, with a propagative part estimated here to be 
18% of the total measured amplitude (based on a gradi- 
ent of magnetic field of 18% along the swimmer's body). 
This model is compared to the experimental results on 
Fig. [5Kb) and shows a reasonable qualitative agreement, 
given the fact that there are no adjustable parameters. 
The main features of the swimming dynamics can there- 
fore be explained by a simple kinematic model based on 
viscous friction associated with the bending modes of an 
elastic rod. Considering the mode 3 above 15 Hz (dash 
line) gives a better agreement than considering the mode 
2 for all frequencies (continuous line); however, the model 
under-estimates the swimming velocity at high frequency. 
An alternative model has been developed, where viscous 
friction is considered for the longitudinal motion and in- 
ertial friction for the transverse motion, without great 
influence on the quantitative results. 

In conclusion, we have designed a novel magnetic swim- 
mer that moves through an original mechanism: genera- 
tion of a bending wave via a twisting/buckling instability, 
and propagation of this wave by breaking the head-tail 
symmetry of the applied twisting torque. The swimming 
dynamics is then mainly governed by the bending modes 
of the equivalent elastic rod. 
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